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Equilibrium atomic configurations and electron energy structure of ethanol adsorbed on the 
Si (111) surface are studied by the first-principles density functional theory. Geometry optimization 
is performed by the total energy minimization method. Several equilibrium atomic configurations 
of ethanol, both undissociated and dissociated, on the Si (111) surface are found. Reaction path- 
ways and predicted transition states are discussed in comparison with available experimental data 
in terms of the feasibility of the reactions occurring. Analysis of atom and orbital resolved projected 
density of states indicate substantial modifications of the Si surface valence and conduction bands 
due to the adsorption of ethanol affecting the electrical properties of the surface. 



I. INTRODUCTION 

Physics of organic molecules adsorbed on semiconduc- 
tor surfaces opens up a new field of applications in high- 
tech industry such as chemical sensors and molecular 
electronic devices [U [2] A detailed understanding 
of the adsorption mechanisms of hydro-carbons on solid 
crystalline surfaces is one of the most important topics 
of modern material and surface science (see e.g. [4 and 
references therein). State-of-the-art first principle meth- 
ods based on the Density Functional Theory (DFT) are 
very useful tools for providing a detailed understanding 
of structural and electronic processes on solid surfaces 

Silicon is one of the basic materials used in mod- 
ern electronics. The atomic processes on Si surfaces 
have attracted interest from both fundamental and ap- 
plied aspects of the surface science Adsorption of 
organic molecules on Si surfaces has been extensively 
studied both experimentally [ll [H [9] and theoretically 
m Uni [11] [12 . Equilibrium atomic structures, thermody- 
namic stabilities, and transition states have been studied 
within last decade [11 [Ml [151 [161 (HI [H |19]. Fourier 
Transform Infrared (FTIR) spectroscopy and Scanning 
Electron Microscopy (SEM) studies demonstrated dis- 
sociative character of organic molecular adsorption on 
Si-surface Due to its applications in micro- 

electronics a great deal of attention was given to the 
study adsorption mechanisms on the Si (001) surface 
[2011211 [221 [231 [24]. 

Ethanol is a relatively small organic molecule, which 
represents both a model system for the investigation of 
organic molecules with silicon surface, but which is also 
of fundamental importance since ethanol is used in many 
processing (cleaning) steps in the preparation of silicon 
surfaces for a variety of technical applications. Nonethe- 
less, little theoretical work has been performed to un- 
derpin experimental studies of the interaction of ethanol 
with Si. Previous ethanol adsorption studies by Sil- 
vestrelli [6j, Eng [9 , and Zhang [13 were done on the 
Si (100) surface concluding that the dissociation kinet- 
ics are energetically favored by 0-H bond cleavage above 
the Si-Si dimer. This adsorption mechanism cannot be 



applied to the Si (111) surface because it is completely 
different from the Si (001) surface with respect to the ge- 
ometric and electronic structure. The Si dimers which in- 
fluence and, to a large extent, define the adsorption mech- 
anisms observed on the Si (001) surface are not present 
on the Si (111) surface. However, the charge separation 
that occurs on the Si-Si dimers, which is also responsible 
for the buckling observed, occurs to and extent on the 
Si (111) surface. The Si-Si dimers, consisting of an up 
nucleophillic and a down electrophillic Si atoms, attract 
through Coulombic forces the ethanol molecule, which 
has significant charge separation of its own, forming a 
dative bond between the down Si atom and ethanol's 
oxygen. At this stage the molecule is providing both elec- 
trons required for the bond. In the presence of ethanol 
similar buckling and charge separation is observed on the 
Si (111) surface, which unlike the Si (100) surface does 
not have charge separation nor buckling in the absence or 
the organic molecule. Due to the Si surface atoms long- 
ing for a stable configuration filling their outer p-shell 
and due to the potential exerted by the ethanol molecule, 
charge separation occurs on the Si (111) surface which 
in turn causes buckling forming a similar configuration 
to the Si (100) surface with a down Si atom immedi- 
ately below ethanol's O-H group and alternating in the 
four directions throughout the surface. This can only be 
observed on an unpassivated Si (111) surface where the 
terminating hydrogens were desorbed previously. 

Nanostructured materials, for example porous silicon 
(PorSi), have experienced a recent surge in interest due 
to their unique properties and applications. PorSi, con- 
sisting of a network of pores of various sizes embedded in 
silicon layer, was discovered to have strong photolumines- 
ence (PL) at room temperature in the visible spectrum 
^25j due to quantum confinement [26 . The luminescence 
from PorSi has been observed to be up to five orders of 
magnitude stronger than that of bulk Si [2T. The en- 
hanced PL emission intensity of PorSi is in part caused 
by surface chemistry [28 which shows strong sensitiv- 
ity to the environmental presence of various molecules 
leading to variations of the PL intensity and peak shift 
[29[ l3Q[ l3l] [32] . impedance [33], capacitance [34., I- 
V characteristics [35^, and Fabry-Perot fringe patterns 
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[361 137] . The large surface to volume ratio of PorSi causes 
the changes on the surface to affect the properties of the 
material as a whole. Thus detailed study of adsorption 
mechanisms on the surface is even more crucial for nanos- 
tructured materials. 

The work presented here is organized as follows: sec- 
tion |TI| describes the computational method. Section 
|III| presents the result separated into three subsections. 
The predicted equilibrium atomic configurations are de- 
scribed in section [Til A[ Transition states along with reac- 
tion coordinate graphs for most favorable reaction paths 
are presented and discussed in section |IIIB| In section 
|III C| the spectra of atom and orbital projected density 
of electronic states are analyzed. 



II. METHOD 

Density functional theory (DFT) has been shown for 
decades to be very successful in the ground state analysis 
of different materials ^j^. For atomic system which in- 
clude molecules and solids it is important to realistically 
describe both short (covalent) and long range (Coulomb 
and van der Waals, vdW) components of intermolecu- 
lar and molecule-surface interactions [39]. Local den- 
sity (LDA) [40] and generalized gradient approximations 
(GGA) [4T are frequently used methods to account for 
exchange and correlation (XC) interaction. The vdW 
interactions is not included in standard DFT. However 
detailed analysis of organic molecule adsorption on solid 
surface demonstrates that around equilibrium the kinetic 
energy of valence electrons remains mainly repulsive, and 
XC effects are mostly responsible for the attraction [39\ . 
It has been shown that equilibrium distance between or- 
ganic molecule and solid (graphene) surface predicted by 
LDA is in good agreement to the value followed from ex- 
plicit inclusion of the vdW into interaction Hamiltonian 

m- 

In this work we used ab initio pseudopotential method 
within the DFT and the supercell formalism to study 
equilibrium geometry, energy kinetics parameters of dis- 
sociative chemical reactions, and electron energy struc- 
ture of Ethanol adsorbed on Si (111) surface. The un- 
reconstructed Si (111) surface is unstable which results 
in its further reconstruction [4 . This surface is also 
rather reactive and hydrogen adsorption is frequently 
used to passivate it. Our adsorption study of ethanol 
on such Si (111) surfaces could be separated into two 
phases: first reactivation of the surface and then the re- 
action of ethanol with dangling bonds. Here we con- 
centrate on the second phase and consider unrecon- 
structed Si (111) surface. The equilibrium geometries of 
CH3CH2OH molecule adsorbed on the Si (111) unrecon- 
structed surface are calculated by LDA method [40 em- 
ploying norm-conserving pseudopotentials [42 using the 
DMol^ [43 software package. We have explicitly checked 
that the structural and binding properties of our system 
are well converged for the double numerical plus polar- 




FIG. 1: (Color) Isolated ethanol molecule. The white, 
red, and grey, spheres represent hydrogen, oxygen, and 
carbon atoms, respectively. 



TABLE I: Predicted structural parameter of the free 
ethanol molecule in comparison with experimental data 



(see Fig. [T]) 




Calculations 


Experiment 1481 


d{0-R) (A) 


0.97 


0.97 


d{c-0) (A) 


1.41 


1.43 


d{C-C) (A) 


1.50 


1.51 


ZC-C-O (deg) 


107.9 


107.8 


ZC-O-H (deg) 


108.5 


105.0 


ZC-C-H (deg) 


109.9 


110.0 


ZC-C-H (deg) 


110.9 


111.0 



C-H bonds have lengths in the range 1.101-1.11. 



ization basis set used. The use of the exact DFT spheri- 
cal atomic orbitals has several advantages. For one, the 
molecule can be dissociated exactly to its constituent 
atoms (within the DFT context). Because of the local 
character of these orbitals, basis set superposition effects 
[44] are minimized and good accuracy is achieved even for 
weak bonds. Absolute binding energy values of the equi- 
librium structures are computed within GGA employing 
ultrasoft pseudopotentials [45 and by using the CASTEP 
[46] software package. The transition states, energy bar- 
riers, and reaction coordinates, have been obtained using 
the "generalized synchronous transit method" [47j. 

For free-standing ethanol molecule the method chosen 
in this work shows excellent agreement with experimental 
results. The structure of the isolated molecule is shown 
in Figjl] while geometrical parameters are reported in 
Table Hr 

The solid surface is modeled as a (2 x 2) periodic struc- 
ture of six Si bi-layers and a 12 A thick vacuum layer 
separating the slabs to eliminate the spurious interac- 
tions between the ethanol molecule and back face of the 
slab; the bottom four Si bi-layers are constrained in or- 
der to reproduce the bulk electronic structure of Si. For 
the surface physics contribution of intermolecular inter- 
action on the surface is important. However, by chosing 
the (2 X 2) unit cell (where adsorbed molecules are far 
apart from each other) this issue remains out of scope 
in this paper. Hydrogen atoms are used to saturate the 
dangling bonds on the back surface of the slab which are 
also allowed to relad providing a measure of interaction 
with the ethanol molecule from the unitcell below. The 
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adsorption energy per molecule is defined as the differ- 
ence between the total energy of the adsorption system 
and the energies of the isolated components, namely the 
clean substrate and the adsorbate, divided by the number 
of adsorbed molecules N per unit cell: 

^ads — {^substrate H~ ^adsorbate 

- Etat)/N (1) 

Note that the adsorption energy is positive by definition 
when the adsorbates bind to the substrate. 



III. RESULTS AND DISCUSSION 

In this section we present first the equilibrium geome- 
tries resulted by atomic relaxation during the adsorption 
of ethanol on Si (111) surface. The second part of this 
section discusses the predicted transition states while the 
related transformations of electronic structure are dis- 
cussed in the third part of this section. 



A. Equilibrium Configurations 

Predicted adsorption energies of Ethanol adsorbed on 
Si (111) are presented in Table [ll| Undissociated and dis- 
sociated atomic configurations are denoted with U and D, 
respectively. Physisorption of Ethanol with CH3 group 
oriented downward to the surface is unfavorable (see Ul 
Fig. [2] and Table [ll| . Consequently the ethanol molecule 
favors the interaction with the Si (111) surface through 
the formation of a dative bond between the O atom and 
a Si surface atom. In these states the ethanol molecule 
remains undissociated with different stable orientations 
with respect to the Si surface (see U2 and U3 in Fig. 
Creation of dative bond results in geometry changes (as 
shown in Fig.[2|, charge redistribution between the O and 
Si surface atoms, and lowering of the total energy (see Ta- 
ble |llj). This reaction does not have an energy barrier and 
occurs spontaneously at room temperature. The C-C/C- 
H bond cleavage has an energy barrier much higher than 
that related to the molecule rotation above the Si (111) 
surface; the ethanol molecule will simply rotate into a 
more favorable, the 'OH group down', configuration. 

The atomic configurations most energetically favorable 
appeared after ethanol dissociation on the surface (see 
Table |llj). It is instructive to consider the energy of tran- 
sition states on the path from undissociated to dissoci- 
ated configurations. Starting from Ul, U2, or U3 states, 
there is a number of possible pathways for continuous 
ethanol reactions with the Si surface, which imply break- 
ing one or more molecule bonds in order to realize more 
stable, dissociative, configurations. We found nine stable 
dissociative configurations which are shown in Fig. [3j 

The Dl and D2 structures are obtained through break- 
ing of the 0-H bond while the detached proton forms a 
new bond with a Si surface atom. Both of these reactions 



have very small energy barriers and are "self-sustaining" 
by drawing energy from the initial adsorption to U2/U3. 
This will be discussed in more details in a subsequent 
paragraph. The D3 configuration is achieved through 
the C-H bond cleavage and forming bonds between Si- 
C and Si-H. The starting point for this reaction is the 
Ul configuration which is energetically significantly less 
favorable than U2/U3. In addition, D3 itself is energet- 
ically less favorable than the previously mentioned Dl 
and D2 configurations. Thus even though the D3 con- 
figuration is possible, it is much more likely to find the 
Dl and D2 configurations. The fourth dissociative con- 
figuration, D4, can be obtained by breaking the covalent 
C-0 bond; the ethyl group attaches to another Si surface 
atom. This is the energetically most favorable configura- 
tion out of all the found dissociated structures which have 
two bonds between surface and molecule. The D5 and 
D6 structures are obtained through cleavage of the C-C 
bond, followed by the methyl and CII2OII groups forming 
new bonds with two neighboring Si surface atoms. The 
starting point for this reaction is Ul. Considering the 
energy required to break the C-C bond and the energy 
of the Ul and D5/D6 structures as compared to other 
configurations, it is highly unlikely that this configura- 
tion will be seen. The D7 structure is obtained through 
the C-C bond cleavage, just like D5 and D6, however this 
configuration follows from the dative undissociative con- 
figurations U2 and U3. The methyl group bonds to a Si 
surface atom, while the proton migrates from the oxygen 
to the danging bond on the carbon. This also converts 
the Si-0 bond to a covalent bond. The two "bridging" 
structures D8 and D9 follow from the D2 and D3 con- 
figurations, cleaving the 0-H bond; two new bonds are 
created between the oxygen and another Si surface atom 
and between the detached proton and yet another Si sur- 
face atom. 



B. Transition States 

The values of activation energy barriers are important 
energy parameters for the surface reaction kinetics. In 
Fig. [4j the transition energy paths are given for a few 
types of the surface reactions. The transition state TSl 
of the reaction U1-D3 is predicted to be above free-energy 
level, the energy of the system where the molecule is in- 
finitely separated from the surface, which indicates that 
this kind of reaction leads to the desorption and it is less 
probable on the Si(lll) surface. 

Note that significant Si host atom reconstruction oc- 
curs in the Ul structure. Instead of the Si surface hav- 
ing bulk-like six- fold rings, the surface has alternating of 
seven- and five-fold rings (see Ul in Fig. [2|. Under ex- 
citation the molecule undergoes a reaction and it arrives 
at D3 product configuration. The D3 configuration is 
more stable than Ul by about 1.4eV. However the en- 
ergy barrier of over 1.42 eV is significant and will not 
occur during normal circumstances. 
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FIG. 2: (Color) Stable configurations for physisorption of ethanol on Si (111). As before, the white, red, and grey, 
spheres represent hydrogen, oxygen, and carbon atoms, respectively, while the yellow spheres represent silicon atoms. 



TABLE II: Binding energies and structural parameters of the configurations shown in Figs, pi and 





Ul 


Undissociated 
U2 


U3 


Dl 


D2 


1 Step Dissociated 
D3 D4 


D5 


D6 


2 Step Dissociated 
D7 D8 D9 


Eads (eV) 


0.93 


1.524 


1.522 


3.25 


3.26 


2.33 


4.00 


2.46 


2.59 


3.57 


4.29 4.25 


d(O-H) (A) 


0.97 


1.05 


1.04 






0.97 


0.98 


0.97 


0.99 






d(C-0) (A) 


1.41 


1.46 


1.46 


1.42 


1.42 


1.41 




1.41 


1.40 


1.41 


1.44 1.45 


d{C-C) (A) 


1.50 


1.49 


1.49 


1.50 


1.50 


1.50 


1.51 










d{C-Si) (A) 












1.88 


1.89 


1.90 


1.93 


1.88 


1.98 1.97 


d{o-Si) (A) 




1.85 


1.85 


1.66 


1.66 




1.67 






1.67 


1.69 1.70 


d(R-Si) (A) 








1.50 


1.51 


1.50 










1.51 1.51 



The activation energy required to go from U2 to D4 
is only leV (see TS2 in Fig. [4|, which is significantly 
less than that required to go for Ul to D3 and it is be- 
low the free-standing energy level and thus can occur at 
room temperature. In addition, the starting structure 
U2 is energetically much more favorable than Ul. The 
resulting configuration D4 is one of the most stable and 
energetically most favorable structures we found. 

The undissociated configuration U2 is transformed 
through dissociation of the Ethanol molecule into D2 con- 
figuration over a very small activation barrier. The en- 
ergy reduction of the transition from U2 to D2 is 1.74 
eV (see Fig. [4|. The predicted energy of the transition 
state TS2 between U2 and D2 is 0.07 eV above U2. This 
indicates that dissociation of hydrogen and its reaction 
with the Si dangling bond is the most probable process at 
room temperature. The molecule shown in D2 can then 
proceed to dissociate further into D8. Though the result- 
ing structure is the most stable structure we found the 
very high energy of the D2^D8 transition step (3.36eV) 
would make this step very unlikely, especially since the 
activation energy for the backreaction D2^U2 is smaller. 
Therefore it would seem likely that dissipation of the en- 
ergy released in the U2^D2 step will cause D2 to become 
a trap state and that even though D8 may have the low- 
est energy state, it will not be reached without significant 
thermal or photo excitation. 



In order to understand the nature of surface electronic 
structure modifications through ethanol adsorption it 
is important to analyze the projected density of states 
(PDOS). Details of the PDOS spectra calculations are 
given in [49]. Electron energy structure changes during 
the U2-D2 reactio n the s- and p-orbital PDOS spectra 
are given in Fig. [5] |(b)| and (c) for both the U2 and D2 
configurations, respectively. 



C. Electronic Structure 

For bulk Si, our GGA-based computational method 
underestimates the gap by about 0.5 eV. A comparison 
of calculated PDOS spectra of U2 and D2 structures with 
those reported earlier for bare and hydrogenated Si (111) 
surface [4 indicates that electronic states in the gap re- 
gion around 0.2 to 0.5 eV correspond to unsaturated dan- 
gling bonds of surface Si atoms. These states are com- 
pletely eliminated after saturation with hydrogen of the 
remaining Si dangling bonds on the surface. As expected 
the saturation of Si-dangling bonds with ethanolic hydro- 
gen due to U2-D2 reaction (see Figs. [2]and|3| results in 
change to the surface energy structure in the gap. Elec- 
tron states changes in the bandgap (see total panel in 
Fig. [5][(b)] and ^ appear due to hybridization caused 
by unpaired Si electrons on the surface (dangling bonds). 




Further, it is important to understand the hybridiza- 
tion changes of electronic bonds on the surface in the 
course of the surface reaction. Atom-resolved PDOS 
spectra provide with additional information about sin- 
gle unit cell atom contributions to the total density of 
states. In Fig. |5] we present both atom, as well as s- 
and p-orbital resolved PDOS spectra related for the re- 
action U2-D2. The U2 to D2 reaction path is energet- 
ically more preferable then other reactions studied here 
(see Table |llj). It is interesting to understand how elec- 
tronic structure is modified by the transformation of the 
bond between the surface and the molecule from dative 
to covalent through dissociation. 

The shape of atom resolved projected density of states 
spectra (AR-PDOS) shown in Fig. [s] reveals that there 
are no electronic states around zero energy related to 
the two ethanol carbons, oxygen, nor any of the carbon- 
bonded hydrogens. The electronic structure correspond- 
ing to these states is located far away from Fermi energy 



and is not discussed here. Fig.[5]shows the AR-PDOS for 
the U2 (dashed black lines) and D2 (solid red lines); only 
atoms that have electron states around eV are shown. 
A comparison of the upper panels in Fig. [5] indicates that 
due to the dissociation of the hydrogen atom and the cre- 
ation of the Si-H bond, the hydrogen related electronic 
states are pushed out of the gap. Unsaturated dangling 
bond states still present in the gap area around 0.4 eV do 
not show rehybridization during the U2-D2 reaction. The 
data in Fig. |5] clearly indicate substantial changes of the 
valence band nature during the U2-D2 reaction: the top 
of the i'-band in D2 configuration is determined mostly 
by the contributions of up and down Si atoms (see atoms 
4 and 5 in Fig. [5] [(b) and (c)), in contrast to a single Si 
atom (hydrogenated Si atom in D2 configuration. Atom 
4) contribution in U2. A similar reconstruction happens 
to the bottom of the c-band (see c-band for atoms ^, ^, 
and 5 in Fig. |5|. The transformation of the Si-0 bond 
from dative to covalent pushes the oxygen related states 



Transition States 



TS4 




FIG. 4: (Color) Reaction coordinates and transition states for typical configurations. As before, the white, red, grey, 
and yellow spheres represent hydrogen, oxygen, carbon, and silicon atoms, respectively. 



out of the gap area (atom 2 in Fig. 5| |(c)[ ) 



For the purposes of photosensitivity, it is instructive 
to consider further desorptive reactions of the ethanol 
on Si-surface. Reaction of U2-D4 results in more ener- 
getically favorable configuration. However the activation 
energy for the U2-D4 reaction (TS2=1.1 eV) is much 
bigger then that for the U2-D2 reaction (see Fig. [4|. 
The predicted configuration of D8 which is characterized 
by almost complete dissociation of the ethanol molecule 
through braking very strong C-C bonds has the lowest 
free energy (see Table [ll[) . The D8 configuration is char- 
acterized by the "bridge" structure, achieved by break- 
ing the C-C and 0-H bonds of Ethanol during dissocia- 
tive chemical reaction on Si(lll) surface. The resulting 
atomic structure bridges two neighboring surface silicon 
atoms leading to the most stable configuration predicted 
in this work. However, this reaction has very high ac- 
tivation barrier (3.36 eV). This energy barrier prevents 
the "bridge" structure from being realized under nor- 
mal conditions. This study clearly indicates that the D8 
configuration-related process could only be possible un- 
der strong (e.^. under laser power) excitations of the 



surface. 

Thus the analysis of our results indicates that the dis- 
sociation of Ethanol on Si(lll) surface, that is accompa- 
nied by the creation of new covalent bonds, is character- 
ized by substantial modifications of the valence band top 
(see the PDOS spectra) . Since top valence band electron 
transitions dominate in the optical absorption threshold, 
one can expect substantial modification of the photosen- 
sitivity of the Si-surface due to the ethanol adsorption. 
The predicted energy characteristics of the Ethanol ad- 
sorption on Si(lll) surface could be used for photoin- 
duced reactions, in particularly of the systems based on 
patterned Si-surfaces (e.g. PorSi) [50^ . 



IV. CONCLUSIONS 

The equilibrium geometries of Ethanol adsorbed on 
bare (initially unreconstructed) Si (111) surface are pre- 
dicted by ah initio pseudopotential method. Desorp- 
tion of ethanol hydrogen and its reaction with the Si- 
dangling bond is shown to be a most probable product of 
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(a) U2 and D2 configurations with individual atoms labeled. 
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FIG. 5: (Color) Atom a nd s hell resolved projected density of states for U2 and D2 con figur ations. Nu mbe rs on 
atomic geometry pictures | (a) | denote surface atoms which correspond to the panels of s- |(b)| and p-shell [(c)] PDOS 
spectra. As before, the white, red, grey, and yellow spheres repr esen t hy drogen, oxygen, carbon, and silicon atoms, 

respectively. The dashed black and solid red lines (figures [(b)] and [(c)] ) represent the U2 and D2 configurations, 

respectively. 



the Ethanol reaction on Si(lll) surface, consistent with 
the results previously obtained for Ethanol adsorption 
on Si (100) surface. More energetically favorable config- 
urations could only be achieved by overcoming substan- 
tial activation barriers. One of the most stable equilib- 
rium structures found is the so called "bridge" structure, 
achieved by breaking two bonds, C-C and 0-H of Ethanol 
during dissociative chemical reaction on Si(lll) surface. 
The resulting atomic structure bridges two neighboring 



silicon surface atoms leading to the most stable configu- 
ration predicted in this work, however it is unlikely that 
this structure can be observed without substantial free 
energy being available in the system causing the D2 to 
become a trap state. The calculated energy barrier for 
the 0-H bond clevage is significantly lower than that ob- 
tained for the Si (100) surface. The absence of the Si-Si 
dimer, which is buckled during the undissociated adsorp- 
tion, benefits the reaction by lowering the energy of the 
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